Recovery of free energy branches in single molecule experiments 
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We present a method for determining the free energy of coexisting states from irreversible work 
measurements. Our approach is based on a fluctuation relation that is valid for dissipative transfor- 
mations in partially equilibrated systems. To illustrate the validity and usefulness of the approach, 
we use optical tweezers to determine the free energy branches of the native and unfolded states of 
a two-state molecule as a function of the pulling control parameter. We determine, within 0.6 feT 
accuracy, the transition point where the free energies of the native and the unfolded states are equal. 
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Recent developments in statistical physics [lj have pro- 
vided new methods to extract equilibrium free energy 
differences in small systems from measurements of the 
mechanical work in irreversible processes (see [51 [3] for 
reviews). In this regard, fluctuation relations |3J are 
generic identities that establish symmetry properties for 
the probability of exchanging a given amount of en- 
ergy between the system and its environment along ir- 
reversible processes. If a system, initially in thermody- 
namic equilibrium, is strongly perturbed by fast varying 
a control parameter A between two values Ao and Ai, then 
the system is driven out of equilibrium. The work exerted 
upon the system, averaged over the ensemble of all pos- 
sible trajectories, reads (W) — (f x * (dH/8X)dX) where 
7i is the system Hamiltonian. According to the second 
law of thermodynamics, (W) is always greater than the 
free energy difference between the initial and final states, 
AG = G(Ai) — G(Ao). The Crooks fluctuation relation 
[I] extends the predictive power of the Second Law by 
establishing a symmetry relation for arbitrary function- 
al of a trajectory T measured along a noncquilibrium 
process (forward or F process) and its time reversed one 
(reverse or R process). In the forward process the system 
starts in equilibrium at Ao and A is varied from Ao to Ai 
for a time tf according to an arbitrary protocol \(t) (i.e., 
Ai = X(tf)). In the reverse process the system starts in 
equilibrium at Ai and A is varied from Ai to Ao following 
the time reversed scheme, given by X(tf — t). In its most 
general form the Crooks fluctuation relation reads [4 

(FeM-(3(W-AG)}) F = (f) R , (1) 

where T stands for an arbitrary functional of the forward 
trajectories the system can take through phase space, /3 
is the inverse of the thermal energy k^T where fee is the 
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Boltzmann constant and T the temperature of the envi- 
ronment. In this relation, T is the time reversal of T, 
while the averages (•)f(r) are taken over the ensemble of 
all possible forward (reverse) trajectories. The particular 
case T = S(W — W(T)) yields a relation between work 
distributions along the forward and reverse processes, 
P F (W) = Pr(-W) exp\J3(W - AG)]. This relation has 
been experimentally tested and used to extract free en- 
ergy differences in single molecule experiments [51|5J[7]. A 
thorough discussion on its validity domain can be found 
in g] 

Fluctuation relation under partial equilibrium con- 
ditions. By considering only the trajectories that go 
from one specific subset of configurations to another one, 
Maragakis et al. [H] have derived another relation use- 
ful to extract free energy differences between subsets of 
states. In principle, the validity of Eq. is restricted to 
initial conditions that are Gibbsian over the whole phase 
space S (what we might call global thermodynamic equi- 
librium). It is, however, possible to extend Eq. (JT|) to 
the case where the initial state is Gibbsian but restricted 
over a subset of configurations (what we might call partial 
thermodynamic equilibrium) . A relation mathematically 
similar to Eq. Q can be derived, but involving nonequi- 
librium processes that are in partial (rather than global) 
equilibrium. It is useful to rephrase here the derivation in 
such a way to emphasize the role played by partial equi- 
librium. As we will see this makes it possible to exper- 
imentally determine the free energy of coexisting states 
for values of A such that the system is never globally 
equilibrated. 

Let P^ q s ,(C) denote the partially equilibrated (i.e., 
Boltzmann-Gibbs) distribution for a given value of A. 
Such distribution is restricted over a subset S' of con- 
figurations C contained in S (i.e., C e S' C S). 
The case S' = S corresponds to global equilibrium: 
P\^(C) = P\s(ty- Partially equilibrated states sat- 
isfy P A cq 5 ,(C) = P a EQ (C)xs'(C)2a,s/Za,ss where X S> is 
the characteristic function defined over the subset S' 
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(xs'(C) = 1 if C £ S' and zero otherwise), and Z\.s> 
is the partition function restricted to the subset S' , i.e., 
z \,s> = Ece5' ex P(~/ 3£: A(C)), with E X (C) the energy 
function of the system for a given A and C. Given a 
forward trajectory T, going from configuration Co when 
A = A to C\ for A = Ai, let So (Si) be the subset of 
iS over which the system is partially equilibrated at Ao 
(Ai). Consider now the following transformation of the 
functional T in Eq. @: T(Y) -> X sM^)Xs x (C r ). 
Under previous conditions the following identity can be 
proved (see Supp. Mat.): 



CFexp[-/?(W-AGf$)])j 



So*— Si 



(2) 



where the average (•)p° ] ^' Sl ^' So ^' Sl ' ) is now restricted to 
forward (reverse) trajectories that start in partially equi- 
librated state Sq (Si) at Ao (Ai) and end in Si (So) at 
Ai (A ). pf'^ Sl (Pr 0< ~ 51 ) stands for the probability to 
be in Si (So) at the end of the forward (reverse) pro- 
cess defined above, and AG^'^ — Gs ± (\i) — Gs a (^o) 
is the free energy difference between partially equili- 
brated states So and Si. In the following, we will drop 
the subscript (F, R), leaving the direction of the arrow 
to distinguish forward from reverse. Moreover, we will 
adopt the shorthand notation Vp = p s °^ Sl / p s °^ s i . If 
T = 1 we obtain a generalization of the Jarzynski equal- 



ity V^(e X p[~P(W - AG£;£)]) = 1. Whereas for the 
particular case T — 5(W — W{Y)), we get the relation 
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exp\/3(W (3) 



l(-W) 



Sq,XqJ 



which has been used in [H] in the case of global equilib- 
rium initial conditions. 

Experimental test. Here we test the validity of Eq. ([3| 
by performing single molecule experiments using optical 
tweezers. Let us consider an experiment where force is 
applied to the ends of a DNA hairpin that unfolds / refolds 
in a two-state manner. The conformation of the hairpin 
can be characterized by two states, the unfolded state 
(U) and the native or folded state (N) — see Fig. [I] The 
thermodynamic state of the molecule can be controlled 
by moving the position of the optical trap relative to 
a pipette (Fig. [2j upper panel) . The relative position of 
the trap along the x-axis defines the control parameter in 
our experiments, A = x. Depending on the value of x the 
molecule switches between the two states according to a 
rate that is a function of the instantaneous force applied 
to the molecule [3]. In a ramping protocol the value of 
x is changed at constant pulling speed from an initial 
value xo (where the molecule is always folded) to a final 
value xi (where the molecule is always unfolded) and the 
force / (measured by the optical trap) versus distance x 
curves recorded. By computing i) the fraction of forward 
trajectories (i.e., increasing x) that go from N (= So) 
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FIG. 1: A two-state molecule in a pulling experiment, (a) 
Schematic picture of the free energy landscape in a two-state 
folder as a function of the reaction coordinate. The blue and 
red colors represent the subsets of configurations that define 
the JV and U states respectively, (b) Depending on the value 
of the control parameter x, the shape of the free energy land- 
scape is tilted toward one state (either N or U). For each 
value of x, the Boltzmann-Gibbs equilibrium value of the 
force restricted to each state defines the native (blue) and 
unfolded (red) force-distance branches, (c) During a ramping 
protocol, x is changed at a constant pulling speed from xq 
to x\ (from x\ to xo in the reverse case) and the molecule 
can visit both states as indicated by the colored circles. In 
order to measure Gu(x) (all free energies are computed with 
respect to Gjv(a;o)) at a given value of x (gray area), only 
forward and reverse trajectories that are in JV at xo and in U 
at x have to be considered (which is true only for the reverse 
trajectory in the example shown in the picture). 



at xo to U (= Si) ai x and ii) the fraction of reverse 
trajectories (decreasing extension) that go from U at x 
to N at xo, we can determine "P^. Then by measuring the 
corresponding work values for each of these trajectories, 
we can use Eq. (|3| to estimate the free energy of the 
unfolded branch Gu(x) as a function of x. By repeating 
the same operation with TV instead of U, the free energy 
of the folded branch Gn(x) can be measured as well. 
Note that we adopt the convention of measuring all free 
energies with respect to the free energy Gn(xq) of the 
native state at xq. We are also able to compute the free 
energy difference between the two branches, AG^(x) — 
Gu(x)-G N (x). 

We have pulled a 20 bps DNA hairpin using a minia- 
turized dual-beam laser optical tweezers apparatus |10j . 
Molecules have been pulled at two low pulling speeds (40 
and 50 nm/s) and two fast pulling speeds (300 and 400 
nm/s), corresponding to average loading rates ranging 
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FIG. 2: (Upper panel) Experimental setup. (Lower 
panel) Unfolding and refolding work distributions P N ^ U (W), 
P Ni ~ u (-W) at x = x = 61.75 nm measured at two pulling 
speeds: 300 nm/s (red -unfolding- and green -refolding-) and 
40 nm/s (blue -unfolding- and orange -refolding-). In the in- 
set we show two force-distance cycles corresponding to low 
(40 nm/s) and fast (300 nm/s) pulling speeds. The gray area 
indicates the mechanical work (= /_ fdx) exerted along the 
green refolding trajectory. Statistics (number of molecules, 
total number of unfolding/refolding cycles): 40 nm/s (2, 223), 
50 nm/s (2, 183), 300 nm/s (3, 337), 400 nm/s (1, 551). More 
details about the data analysis procedure can be found in the 
Supp. Mat. 



FIG. 3: Experimental verification of the fluctuation re- 
lation Eq. O). (Upper panel) Plot of logP^(^) + 
log{P N ^ u (W)/P N ^ u (-W)) as a function of W (in k B T 
units) at different pulling speeds. A least squares fitting 
method of the experimental data to straight lines gives the 
following slopes: 0.91(4) (400 nm/s), 1.07(8) (300 nm/s), 
1.02(6) (50 nm/s), 1.07(3) (40 nm/s). (Lower panel) Ben- 
nett's acceptance ratio method. We plot the function z(u) 
Eq. Q for different pulling speeds. The thick black line cor- 
responds to the function z — u. Error bars correspond to the 
standard deviation of the Bennett estimate, as determined in 
Ref. [11| . Data statistics are reported in the caption of Fig. [2] 



between 2.6 and 26 pN/s, from xq — to x\ — 110.26 nm 
(for convenience we take the initial value of the relative 
distance trap- pipette equal to 0). A few representative 
force-distance curves are shown in Fig. [2] (inset of lower 
panel). We have then selected a value x = x — 61.75 nm, 
close to the expected coexistence value of x where both 
N and U states have the same free energy (see below). 
Such value of x is chosen in order to have good statistics 
for the evaluation of the unfolding and refolding work 
distributions. The system is out of equilibrium at the 
four pulling speeds. To extract the free energy of the un- 
folded branch Gjj(x), we have measured the work values 

W N^U = J-* f dx ^ W N<-U = f dx abng the unfold 

ing and refolding trajectories, respectively, and then de- 
termined the distributions P N ^ U (W) and P N< ~ U (-W). 
In the main panel of Fig. [2] we show the work distribu- 
tions obtained for a slow and fast pulling process. Note 
that the support of the unfolding work distributions is 
bounded by the maximum amount of work that can be 
exerted on a molecule between xq and x. This bound 
corresponds to the work of those unfolding trajectories 
that have never unfolded before reaching x. 



As a direct test of the validity of Eq. ([3]), in 
the upper panel of Fig. [3] we plot the quantity 
log(P N ^ u (W)/P N ^ u (-W))+logV%(x) against W (in 
k&T units). As expected, all data fall into straight lines 
of slope close to 1. The intersections of such lines with 
the VF-axis provide an estimate of Gjj(x). Note that both 
fast and slow pulling speeds intercept the horizontal axis 
around the same value within 0.75 fee? 1 of error. With 
this method, we estimate Gjj{x) — 185.9(5) k^T. Note 
that this free energy estimate, as all the others in this 
paper, refers to the whole system, comprising hairpin, 
handles, and trap. 

A more accurate test of the validity of Eq. |3| and a 
better estimation [TT] of the free energy Gu(x) can be 
obtained through the Bennett acceptance ratio method 
[L2] . In Bennett's method we define the following func- 
tions: 
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FIG. 4: Free energy of the folded and unfolded branches. The 
purple line is Gn(x), the dark green line is Gu(x), both com- 
puted with respect to Gn(xo). The black solid line and the 
gray area mark the expected value x c — 62.0(7) nm and its 
standard error, while the dashed black line is the crossing 
point x c « 61.28 nm. Both left and right insets show a mag- 
nified image of the coexistence region. In the right one, the 
dashed lines represent the free energies that we would obtain 
if we had overlooked the factor V%. The top panel shows 
the functions AG^i) - AG^(i) and AG^i) - AG"(i), 
which, by definition, should be independent of x and equal to 
Gu(xi). Again, the dashed lines represent the result if we do 
not include the factor V%. 



From Eq. ([3| it has been proved [TTJ [13] that the solution 
of the equation z{u) = u, where 



z(u) = zrefO) - zunf(w) - \ogV%{x) 



(4) 



is the optimal (minimal variance) estimate of Gu(x). 

In Fig. [3] (lower panel) we plot the function z(u) 
for different pulling speeds. It is quite clear that the 
functions z(u) are approximately constant along the u 
axis and cross the line z = u around the same value 
Gu(x) = 186.0(3). A distinctive aspect of Eq. ^ is the 
presence of the factor V^. If such correction was not 
taken into account then the fluctuation relation would 
not be satisfied anymore. We have verified that if Pjy is 
not included in the analysis, then the Bennett acceptance 
ratio method gives free energy estimates that depend on 
the pulling speed (see Supp. Mat.). 

Free energy branches. After having verified that the 
fluctuation relation Eq. ([3| holds and that it can be used 
to extract the free energy of the unfolded (Gu(x)) and 
folded (Gjy(a;), data not shown) branches, we have re- 
peated the same procedure in a wide range of x values. 
The range of values of x is such that at least 8 trajectories 



go through N (or U) (ensuring that we get a reasonable 
statistical significance). The results of the reconstruc- 
tion of the free energy branches are shown in Fig. |1J The 
two free energy branches cross each other at a coexis- 
tence value x c at which the two states (N and U) are 
equally probable. This is denned by AG^(x c ) = 0. We 
get x c ~ 61.28 nm, which is in good agreement with an- 
other estimate, x c — 62.0(7) nm, obtained interpreting 
experimental data according to a simple phenomenolog- 
ical model (see Supp. Mat.). In the insets of Fig. [4] 
we zoom the crossing region. It is interesting to recall 
again the importance of the aforementioned correction 
term (VYj) to Eq. If such term is not included in the 
analysis then the two reconstructed branches never cross 
(bottom right inset). This result is incompatible with 
the existence of the unfolding/refolding transition in the 
hairpin, showing that the factor is key to measure 
free energy branches. 

Equation ([2| is valid in the very general situation of 
partially equilibrated initial states which, however, are 
arbitrarily far from global equilibrium. This makes the 
particular case Eq. (§ a very useful identity to recover 
the free energy of states that cannot be observed in con- 
ditions of thermodynamic global equilibrium. We have 
shown how it is possible to apply Eq. ^ to recover free 
energy differences of thermodynamic branches of folded 
and unfolded states in a two-state DNA hairpin. These 
methods can be further extended to the recovery of free 
energies of non-native states such as misfolded or inter- 
mediates states. 
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I. CROOKS RELATION IN PARTIAL THERMODYNAMIC EQUILIBRIUM 

In this section, we provide a full derivation of the relation (2) that has been introduced in the main article: 



Sa^Si pSo^Si 



(1) 



where the averages (^So^SiiSo^S!) arc taken ovcr a \\ forward (reverse) trajectories that start in partially equilibrated 
subset of state Sq (Si) at A (Ai) and end in Si (So) at Ai (A ). Because initial conditions of the forward and 
reverse processes are taken in partial equilibrium, Eq. (1) allows to compute the free energy differences (AG^J'J° = 
G ! 5 1 (Ai) — Gs (Ao)) between subsets of states in far from global equilibrium conditions (see main text). 

In relation (1), T stands for an arbitrary functional of the forward trajectories the system can take through phase 
space, T is the time reversal of T, W is the mechanical work exerted upon the system along the forward trajectories, 
and p s <>^ s i (p So ^ Sl ) stands for the probability to be in Si (S ) at the end of the forward (reverse) process. In what 
follows, • always refers to the time reversal of •. 

In order to prove the validity of this relation, we adopt a framework similar to the one used for the derivation of the 
generalized Crooks relation [1] (Eq. (1) of the main text). In this regard, we recommend [1] for a thorough discussion 
about the generality of the derivation. 

Thus, we consider a system S composed of a finite set of microscopic configurations C, which evolves according to 
a discrete time Markovian dynamics that verifies the following detailed balance equation: 

p(c! + T^cJ = ex P{^[^(*+D( C '+i) - E Ht+i)(C t )]} (2) 

where p stands for the probability transition between two successive configurations in time, that is C t and C t +i here. 
E\(C) is the energy of S when the control parameter is equal to A and the microscopic configuration is C. This models 
the evolution of a system that can exchange energy with a constant temperature heat bath. We define the forward 
trajectory 7c = {Co, Ci, Ct f -i} so that A(0) = A and X(tf) — Ai. V(Tc) = Ilt^o 1 P(^t — * Ct+i) is the corresponding 
probability for the trajectory 7c to occur during the forward process. Then by using the detailed balance condition, 
we get: 

v \ T c) IltLo P( c t c t+i) llt=o PK c t+i -» C t ) y t=0 J 

where we have used Ct = Ct f -t- At this point, it is useful to remind the definition of the thermodynamic quantities in 
this Markovian process. During the forward process, for a given trajectory 7c one can define the total heat exchanged 
with the thermal bath, Q(Tc), and the total work exerted upon the system, W(T C ), as [2]: 

W(T C ) = J2l E Ht+i)(Ct) - E Mt) (Ct)} 
t=o 
tf-i 

Q(T C ) = £[£ A(t+1) (C t+1 ) - £A(t+i)(C t )] (4) 
t=o 

As a result, eq.(3) can be written as V(Tc)/V(Tg) = exp(— (3Q(Tc)). Moreover by multiplying the equality (3) both 
at the denominator and at the numerator by the respective global equilibrium probability of the initial state, one gets 
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P^ Q (C )V(Tc)/P^(Co)V(f d ) = cxp(f3W{T c ))Z Xu s/Z\ ,s where Z x ,s is the partition function of S. Next, a useful 
way to restrict this relation to partial equilibrium (sec main text) is to introduce the characteristic function defined over 
any subset S' of S: XS'(C) = 1 if C € S' and is equal to zero otherwise, so that P^ q s ,(C) = Pf®(C)xs'(C)Z\,s/Z\,s r 
where P^ q s ,(C) is the Boltzmann weight of C restricted to the subset of state S' G S. Then, from relations (3) and 
(4), we get: 

P \o,S ( C o)' P (' I c) = Z Xl , Sl Xs (Cg) c f 3W (Tc) (5) 

Now, (•) 5 "^ ,s i is an average over the forward trajectories that start in partially equilibrated states So at A and 
end in Si at Ai. Thus, for any functional T of these trajectories, we have: 

E P Zs (Co)V(T c ) . X5l (C tf ) E PZ So (Co)V(T c ) . X5l (C tf ) 
T, P Zs a (CoMT c ) XSl (C tf ) P So ^ 

Tc 

where the sum is over all the possible trajectories of the system. In this relation, we have defined p s °^ s ^ — 
E -P^g n (C )'P(7c)x5 1 {Ctj), that is the probability to be in S\ at the end of the forward process, i.e. starting in 

Tc 

partial equilibrium in So- Applying the relation (6) to jFexp[— /3(W — AG s ^' Xo )] and using relation (5), we successively 
get: 

Y. P tls a {Co)V{T c )Hrc)eM-m - AG*£)] X5l (C t/ ) 
(TeM-m-^Gf±)]) s ^ = * -g-^ (7) 



E P Z Sl {Ci)V{t £ )HTc) exp[/?AG|-^ )] f^XSo (Co) 

Y, p Z.sA)n f c)H%)xs (c tf ) 

% 

E P A7,5 1 (^)^(t-)^)X5 (C; / )E P A7,5 1 (^)^(t-)X5 (C t/ ) 

% % 

T, p Zs^)n f c)xs (Ct,) p So 



(8) 
(9) 



.s^ Sl - (10) 



where from (7) to (8), in addition to the relation (5), we have used XSi(Ct f ) = XSi(Co)- From (8) to (9), we have used 
exp[/3AGg^'^)] Z ^' S g = 1) XSo(Co) = X5 (C* / ), the definition of the time reverse of J 7 , that is T(Tq) = T(Tc), and 
the fact that the set of reverse trajectories is the same as the set of forward trajectories. From (9) to (10), we have 
simply multiplied by the same quantity both the denominator and the numerator. Finally, from (10) to (11) we have 

used the equivalent of relation (6) for the reverse process plus the definition p s i^ s " = ^ P Z Si (^i)^C^j)XSo(^t/) ■ 

% 

This ends the derivation of relation (1). 
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FIG. 1: Plot of log(P N ^ u (W)/P N ~ u (-W)) as a function of W (in k B T units). Compare to the upper panel of Fig. 3 in the 
main text. 
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FIG. 2: Bennett's acceptance ratio method, applied to the erroneous free energy AG, obtained by imposing = 1 in Eq. (12). 
This picture may be compared to the lower panel of Fig. 3 in the main text. 
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II. FURTHER DETAILS ABOUT THE EXPERIMENTAL DATA 

A. Effect of the term V% on the validity of the fluctuation relation 

One of the goals of our research was to provide experimental evidence of the validity of the equation (labelled (3) 
in the main text of the Letter) 

ps^s^w) = -Mm - AGf±)] . (12) 

A consequence of the above relation is that, since the term = p s «^ s ^ /pSo^Si j g obviously influenced by the 
pulling speed, if we ignore it and erroneously apply the usual Crooks relation 



exp[/?(W-AG*£)], (13) 



we are bound to find an incorrect "free energy" difference AG, the value of which depends on the pulling velocity. 
This dependence is illustrated in Fig. 1 and, even more clearly, in Fig. 2. In Fig. 1 the slopes of the fitting lines 
strongly deviate from 1 as the pulling speed changes, while Fig. 2 shows how the experimental lines intersect the line 
z = W at different values that are separated up to 2 k^T. More important, data for fast and low pulling speeds 
(orange and red lines versus blue and red lines, respectively) are clustered in two distinct regions in the figure. This 
shows that the free energy difference AG systematically changes with the pulling speed. 



B. Note on the data analysis procedure 

The work probability density functions (PDF) represented in Fig. 2(b) of the Letter have been calculated by 
following a histogram- free recipe recently proposed by Berg and Harris [3]. The method, in a nutshell, consists of 
approximating the empirical cumulative distribution function (ECDF) with a diffcrcntiablc function that incorporates 
a suitable number of terms of the Fourier expansion of the data. The Kolmogorov test is used to decide when the 
difference between the ECDF and its analytical approximation may be ascribed to chance. The PDF is then evaluated 
simply by differentiation, while error bars can be estimated using a jackknifc method [4]. 

This procedure is especially convenient for our data, because the statistics of the forward trajectories such that the 
hairpin is unfolded at x — x c can get so sparse that the traditional, histogram-based method becomes excessively 
dependent on the number and origin of the bins. A small downside, on the other hand, is that Berg and Harris' 
technique (like the Kolmogorov test it is based upon) is more sensitive in the central region of the distribution than 
about the tails. Since we use the analytic values of the PDF to compute the probability ratios used for Fig. 3(a) in the 
Letter (or Fig. 1 in the appendices), the slight misrepresentation of the extrcma of the unfolding work distributions 
results in the small curvature that can be noticed in those figures. 



III. DERIVATION OF x c FROM A SIMPLIFIED MODEL 

Let us consider a pulling experiment where a two-state-like molecule is driven from an initial configuration charac- 
terized by the control parameter xq <C x c to some final x\ ^> x c , where x c is the value of x for which the states ./V and 
U have same free energy. In this section, a simplified version of the theoretical model described in [5] is introduced, 
that allows the estimation of x c in terms of quantities easily measured in single molecule experiments. 

Let us assume that the Thermodynamic Force-Distance Curve (TFDC, i.e., the force-distance curve that we would 
observe if we were performing the experiment in equilibrium conditions) is shaped like in Fig. 3, that is, it jumps from 
the line f(x) = f Q + k e gx to the line f(x) = f + k e gx — Af at some coexistence point x = x c . This is not the most 
general case, but is a useful approximation for the experimental setting we worked with. 

The coordinates of the points in Fig. 3 are therefore 

A=(X , f + k eS X = /min) , G = (x c ,f + k cS x c ) , 

D={x c , f + k cff x c - Af) , F = (x!,f + fccffxi - A/ = / max ) . 
From the experimental data, it is easy to measure directly fo, k e s, and Af. In the case of our experiment on the 
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FIG. 3: Schematic view of the Thermodynamic Force-Distance Curve. 



two-state hairpin, for instance, we have 

fo = 10.86(1) P N, 
fc off = 48.3(5) pN//im, (14) 
Af = 0.92(3) P N, 

where the uncertainty is estimated as standard error. Moreover, performing many pulling experiments and applying 
the Crooks relation [1] to the resulting forward and reverse work distributions, we are able to determine the reversible 
work W lev , that corresponds to the entire area below the TFDC from x = x to x = X\. From the knowledge of the 
three parameters (14) and of the reversible work W Iev = 353.3(4) k^,T , it is possible to determine x c as we explain 
below. 

First, we introduce the notion of coexistence force / c , defined as the force at the midpoint between C and D. Then, 
we baptize B and E the points where the coexistence force is reached on the folded and the unfolded line, respectively. 
In other words, 

B=(x c - Ax/2, f c ) , E= (x c + Ax/2, f c ) , 

where 

f c = fo + k sS x c -Af/2, A^=^. 

The area (red in Fig. 3) below the TFDC between B and E is what in [5] is called W r i p . The blue area is (fo + 
k e sxo)(xi — xq — Ax), while the green area is k e g(xi — xq — Ax) 2 /2, hence 

Wrip = Wrov - (fo + k c f i X )(x 1 - x - Ax) - fc e ff (^l - x - Ax) 2 /2 . (15) 

On the other hand, it is also true that 

W rip = f c Ax. (16) 
Combining Eqs. (15) and (16), and solving with respect to x c , we find 

W rev - (Xl - X )(f + k eS X ) (x X - Xq - Ax) 2 Ax 



Af 2Ax 2 



+ x , 
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that is (using our experimental conditions T = 24°C, x n = 0, and x\ = 110.26 nm), 

x c = 62.0(7) nm. 

This may be compared to our experimental curves for the free energy, which cross each other around x w 61.28 nm. 
As a consistency check of the model we are using, we can notice that the parameter Ax that represents the length of 
the unfolded hairpin (actually, the difference in length between the folded and unfolded states) is 19.0(6) nm, while 
the coexistence force is 13.4(3) pN, both values being consistent with hopping experiments carried out with the same 
molecule (unpublished data). 
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